clc
h = 0.04;
hold on
while h < 0.26
    c = [0,0,0,0,0,1];
    A = [0,0.025,0,0,0,-1
        0,0,0.015,0,0,-1
        0,0,0,0.055,0,-1
        0,0,0,0,0.026,-1
        -0.05,-0.27,-0.19,-0.185,-0.185,0];
    b = [0,0,0,0,-h];
    Aeq = [1,1.01,1.02,1.045,1.065,0];
    beq = 1;
    LB = zeros(6,1);
    [x,Q] = linprog(c,A,b,Aeq,beq,LB);
    Q = -Q;
    plot(h,Q,'db');
    h = h + 0.01;
end
xlabel('h'),ylabel('Q')
